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Scientific Objectives 


1 . To establish the relative contributions of gravity-driven fluid 
flows, both in liquid and vapor, to ( 1 ) the non-uniform 
incorporation of impurities and defects and (2) the deviation 
from stoichiometry and (3) the compositional variation 
observed in the grown crystals. 

2. To assess the self-induced strain developed during processing 
at elevated temperatures and retained on cooling caused by the 
weight of the crystals. 

3. The relation between fluid phase processes and the generation 
of defects in a grown crystal is an outstanding problem in 
materials growth. Studies in microgravity will be compared 
with modeling and will lead to a greater understanding of the 
processes involved. 


Crystal Growth Activities at MSFC 


Technological significance: 

Growth (melt and vapor) and characterization of II- VI 
compounds semiconducting materials, such as 
HgCdTe, HgZnTe (for IR detectors), CdS and ZnO 
(for UV detector), ZnSe, ZnSeTe (for green /blue 
laser), CdTe and CdZnTe (for x-ray, gamma ray 
detectors). 
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II- VI semiconducting compounds grown at MSFC 


Compounds 

HgTe HgCdTe CdTe CdZnTe ZnTe CdS ZnSe ZnS 

Melting 
points (°C) 

670 700 1092 1130 1292 1397 1526 1718 

Melt growth 



PVT growth 
temperature (°C) 


850 


1000 985 1120 1150 



Crystal growth by Physical vapor Transport 



Results of Cl doped CdTe crystal grown by PVT in microgravity on the 
EUREKA-1 mission. The plots are two dimensional resistivity distributions 
mapped across the axial direction of two crystal wafers. The wafer on the left 
was grown in the microgravity environment, while the wafer on the right was 
grown on Earth. (Benz, et al. 1997) 
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Mass Flux (mole/cm 2 -sec) 


One-Dimensional Diffusion model of PVT 
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a(L)=1.9 



AT (°C) 


Calculated mass flux of 
ZnSe as a function of 
AT and different values 
of a(L). The source 
temperature was 
1080°C. Solid lines are 
for a(L) > 2 and dashed 
lines are for a(L) < 2. 



Calculated mass flux of 
ZnSe as a function of 
AT under the same 
conditions except a 
residual gas pressure 
of 0.008 atm is present 
in the system. 



Calculated mass flux of 
ZnSe as a function of 
residual gas pressure 
for source temperature 
at 1080°C and different 
values of a(L). 


Su, et al. (1998); Sha et al. (1995) 


Summary of one dimensional diffusion analysis 


• Four experimentally adjustable parameters, the source 
temperature, the deposition temperature, the partial pressure 
ratio over source and the residual gas pressure, determine the 
diffusive mass flux of a PVT system. 


• However, two of these four parameters, the partial pressure 
ratio over source and the residual gas pressure, are more 
critical than the others. As will be shown, these two 
parameters are critically dependent on the proper heat 
treatments of the starting materials for optimum mass flux. 


Heat Treatment of CdTe Starting Materials 


Table 2. Amounts of CdTe loaded, heat treatment conditions, the measured partial pressure 
of Te 2 , P Te , and the calculated pressure ratio, p C( /p Te , at 870 ‘ c for as-synthesized measured partial pressure of 

2 2 Te 2 for cell 10A 


and heat treated batches. 


cell 

amount of 
CdTe (g) 

heat-treated 

temperature 

(°C) 

heat-treated 

duration 

(min) 

P T e 2 ( atm ) 
at 870 °C 

p cd/P Te 

1C 2 

at 870 °C 

6 

15.0 

- 

- 

4.06xl0' 3 

1.92 

7 

15.0 

- 

- 

1.65x1 0' 1 

0.00742 

7A 

13.0 

870 

8 

3.44x1 0* 3 

2.47 

9 

12.4 

“ 

- 

8.90xl0' 2 

0.0187 

10 

14.9 

- 

- 

4.95x1 0' 3 

1.43 

10A 

20.3 

870 

8 

2.74x1 0‘ 3 

3.47 

11 

6.9 

“ 

- 

1.51xl0' 2 

0.268 

11A 

20.0 

860 

8 

2.96x1 O' 3 

3.11 

13A 

17.9 

800 

8 

4. 18x1 O' 3 

1.84 

i 



Su et al. (1998) 



Optimum Heat Treatments of Starting Materials 


P Se in equilibrium with ZnSe for various optical cells 

H 2 reduction and vacuum treatments forZnSe(0)-3 established to be the optimum procedure 

1145 °C 1035 °C 



1000/T (K) 


Su et al. (1998) 
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Physical Vapor Transport: 
temperature profile and initial ampoule positions 

• The growth ampoules can be equipped with optical windows to confirm the 
stoichiometry of the starting material before growth. 

• The thermal profile, with a maximum in the middle, was provided by a three-zone 
furnace with an adiabatic zone between central and cold zones. 

• Growth was initiated by translating furnace to the right. Gravity 




Distance 


Su (1995) 
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Steps of PVT Seeded Growth of CdTeS-14 



2. ampoule after soaking at 
saturated position for 12 hours 
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3. reinsert ampoule and process PVT growth 



The initial position of the growth ampoule 
needs to be carefully located so that (1) to 
avoid spontaneous nucleation (2) to 
sublime part of the single crystal seed but 
(3) not completely dissolve the seed. 
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PVT Growth of Various Semiconductors 



seeded growth of CdS 


unseeded growth of PbSe 


unseeded growth of CdS 


unseeded 
growth of ZnTe 
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PVT Growth of ZnSe and Related Ternary Compounds 


1 . Self-seeded growth of ZnSe in vertical (stabilized and 
destabilized) and horizontal configurations 

2. Seeded growth of ZnSe in vertical and horizontal 
configurations 

3. Self-seeded growth of Cr-doped ZnSe in vertical and 
horizontal configurations 

4. Self-seeded growth of ZnSeTe in vertical and horizontal 
configurations 

5. In-situ and real-time optical monitoring of seeded growth 
in a horizontal configuration 
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Flow chart of sample characterization plan 
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Gravity Effects on the Grown Crystals 


Effects were studied by comparing the following 
characteristics of horizontally and vertically grown ZnSe 
crystals in : 


• Grown crystal morphology : contactless growth for the 
horizontal configuration. 

• Surface morphology of the grown crystals was examined by 
SEM and AFM. (growth was terminated by stopping furnace 
translation, lowering the source temperature by 1 0 °C and then 
cooling the thermal profile at the same rate) 

• Segregation and distribution of defects and impurities in the 
grown crystals was determined by photoluminescence, SIMS 
and precision density measurements. 
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Gravity Effects on the Grown Crystals 


Morphology of the as-grown crystals: 

I. Self-seeded ZnSe: Crystals grown in the horizontal configuration grew away 
from the ampoule wall with large (110) facets tend to align parallel to the 
gravitational direction. Crystals grown in the vertical configuration grew in 
contact with the wall to the full diameter. 

II. Seeded ZnSe: the as-grown seeded crystals for the horizontal and vertical 
configurations showed similar characteristics in the morphology as described 
above for the self-seeded growth. 

As-grown surface morphology: 

I. As-grown surface of horizontally grown crystals was dominated by (110) 
terraces and steps (identified to be (221) in one case). 

II. As-grown surface of the vertically grown 

(a) Crystals showed granular structure with nanotubes (200nm OD, 75nm ID, 
25nm in height for one case on ZnSe) on the top. 

(b) Some crystals showed a network of high plateau with each island 30 - 
70mm in diameter and 3.5mm in height. Numerous nuclei were observed with 
diameter 20 - 50nm and height of 1 - 7nm on top of these islands. 
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Morphologies of Self-seeded Vertically 
Grown ZnSe Crystals 




ZnSe-25: vertically stabilized ZnSe-3 1 : vertically destabilized 

ampoule ID: 15mm 

- 1 - 
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Morphologies of Self-seeded Horizontally 

Grown ZnSe Crystals 
ZnSe-44 ZnSe-43 
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Gravity Effects on As-grown crystal morphology 


• For furnace translation rates higher than the mass flux: 

- In the horizontal configuration the crystal maintained the growth surface 
morphological stability by (1) self-adjusting the degree of supersaturation 
to increase the mass flux or/and (2) by reducing the cross section area of 
the grown crystal. 

- In the vertical configuration the crystal growth surface became 
morphologically unstable with voids and pipes embedded. 


ZnSe-47H 7.5 mm/day 


ZnSe-35V 11.4 mm/day 




Morphologies of seeded vertically grown ZnSe 

crystals 



g 


ZnSe(S)-9V 


ZnSe(S)-12V 





Morphologies of the seeded horizontally grown 

ZnSe crystals 


ZnSe(S)-llH 


8 x 


ZnSe(S)-13H 
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ZnSe(S)-8H 



g 




Gravity Effects on Impurities and Defects Distribution 


I. Results from SIMS mappings : 

For the horizontally grown self-seeded ZnSe crystal [Si] and [Fe] showed clear 
segregation toward the bottom on the wafer cut axially along the growth 
axis. 

For the vertically grown seeded ZnSe crystal [Si] and [Cu] showed segregation 
toward the edge of the wafer cut perpendicular to the growth axis . 

II. Mappings of PL near band edge intensity ratios indicated: 

(1) All the horizontally-grown crystals showed the following trends in the 
radial and axial segregation of [Al] and [V Zn ] : 

[Al] segregates radially toward the top and axially toward the first grown 
region and [V Zn ] segregates radially toward the bottom and axially toward the 
last grown region. 

(2) The as-grown surface of the seeded vertically stabilized grown crystal 
showed [Al], [Li or/and Na] and [V Zn ] segregate radially toward the center. 

(3) The as-grown surface of the self-seeded vertically destabilized grown 
crystal showed [Al] and [V Zn ] segregate radially without an apparent pattern. 
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Su et al (1999) 



Impurities Distribution by SIMS (horizontally grown) 



13 - 



11 - 



cone. Si (ppb) 

2.04312 

1.96625 

1.88937 

1.8125 

1.73562 

1.65875 

1.58187 

1.505 

1.42812 

1.35125 

1.27437 

1.1975 

1.12062 

1.04375 

0.966875 


i i i i I 
10 


y (mm) 






Impurities Distribution by SIMS (vertically grown) 


141- 


on of SI In ZnSe- 16 v 



10 11 12 

x (mm) 


oonc. Si (ppb) 

6.66937 

6.23875 

5.80812 

5.3775 

4.94687 

4.51625 

4.08562 

3.655 

3.22437 

2.79375 

2.36312 

1.9325 

1.50187 

1.07125 

0.640625 
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15 


X 

gravity and 
growth direction 


In ZnSe- 16 v 


24 



cone, of Fe (ppb) 

13190.5 
12528 

11865.5 
11203 

10540.5 
9878 

9215.5 
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7228 
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5240.5 
4578 
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Concentration of Cu In ZnSe- 16 v 
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cone, of Cu (ppb) 

1015.56 

987.125 

958.688 

930.25 

901.812 

873.375 

844.938 

816.5 

788.062 

759.625 

731.188 
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645.875 
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Summary on the l 2 and l 1 deep emissions 


l 2 emission: 

• l 2 , the exciton bound to substitutional donor, emission in our ZnSe samples can 
be attributed mainly to Al impurity, with A(I 2 )/A(F X ) = 4.88 corresponding to 1700 
ppb, atomic, or 7.46x1 0 16 cm -3 . 

• Isshiki et al. (1991) gave the expression between intensity ratio (I 2 )/(If x ) and N D : 

log^lAx) = -22.0775 + 1.46268 log 10 N D (cnv 3 ) 

• Therefore, (l 2 /l Fx ) = 82 [A(I 2 )/A(F X )] 


l., dee P emission: 

• I^^p is related to exciton bound to V Zn deep acceptor and [V Zn ] is proportional to 
A(l 1 c,ee P ) / A(F x) . 

• The reaction during Zn vapor annealing: 5000 

Zn zn V z „ + Zn(g) | 

K = [V Zn ] x P Zn = K 1 A(l 1 dee P)/A(F x ) x P Zn J 

£ 2000 

• The ZnSe samples were annealed at 1104 °C: | 

lOOO 

1. A(l 1 dee P)/A(F x ) = 7.52 when sample is in equilibrium with P Zn 

= 6.1x1 0 -3 atm (a = 6.05) A(l 1 dee P)/A(F x ) x P Zn = 0.0459 atm. Photon Energy (eVl 

2. A(l 1 dee P)/A(F x ) = 5.18 when sample is in equilibrium with P Zn 
= 9.0x1 O' 3 atm (a = 19.43) A(l 1 dee P)/A(F x ) x P Zn = 0.0466 atm 
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Distribution of [V Zn ] and [Al] in ZnSe (horizontally 

grown) 



[Al] 


9.5 
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Distribution of [V Zn ] and [Al] in ZnSe (horizontally 

grown) 
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Distribution of [V Zn ] and [Al] in ZnSe (vertically 

grown) 
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Distribution of [V Zn ] and [Al] in ZnSe (vertically 

grown) 


Vertical destabilized configuration 
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Numerical Modeling of 
Physical Vapor Transport 


• Two dimensional and three dimensional calculations 

• Finite element technique - Fidap code 

• Thermal and Species induced buoyancy forces 

• Compressible or Boussinesq model 

• Benchmark -2D (H 2 -l 2 system - PVT growth) 

• ZnSe calculations with residual gas 

• Benchmark -3D (Natural convection in a cylinder) 

• 3-D ZnSe calculations with residual gas 

• Ongoing and future work 
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Ramachandran et al. (2000) 


The Physical Model 


gravity 


Ts 

C 


AS 


source 



Wall eci8f = o 


crystal 


wall 

temperature 


Benchmark 


the H 2 -l 2 System 


• benchmark calculations by Rosenberger et al. (J. Crystal 
Growth 51 426 1981; 67 241 1984; 118 49 1992) 

• source temperature Ts=370.5 K 

• crystal temperature Tc=358.1 K 

• ampoule pressure : 100 torr 

• I . (M=254) is the deposited species and H 2 (M=2.016) is inert 

• 2-d Cartesian system 

• linear wall temperature 

• quasi-compressible and Boussinesq calculations 

• Peclet number analysis, Pe ~ 1 for diffusive flow 
32 • growth rate results 



Parameters of ZnSe System with Residual 

Gas 


a = 2.9 

• density, 1.2xl0' 5 g/cm 3 

• dynamic viscosity, 4.3xl0' 5 Pa-s 

• kinematic viscosity, 36 cm 2 /s 

• diffusivity, Zn in N 2 , 64.59 cm 2 /s 

• diffusivity, Se 2 in N 2 , 71.46 cm 2 /s 

• thermal expansion coefficient, 7.1xl0‘ 4 K' 1 

• Prandtl number, 0.439 

• Schmidt number, Sc zn = 0.557 

• Schmidt number, Sc se = 0.503 
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Velocity Difference Plots - procedure 


• calculations for pure diffusion limited conditions, Og 

• calculations for other g conditions, a conditions, etc. 

• calculate differences in axial (u) and transverse (v) 
velocities at all identical grid locations between previous 
steps; e.g. U(0g) - U(lg) and V(0g)-V(lg) 

• contour the velocity differences and plot 

• determine g-sensitivity 

• Note: all calculations were using the Boussinesq approx. 

I 


conventions: 



y,V f 


source 



crystal 

Tc ; [T=T( Y )] 


x,U 


0=0; horizontal 

0=90; vertical stabilized 

0=270; vertical destabilized 
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ravity Effect on Velocity Difference (horizontal 

case) 


- g-effects with constant crystal temperature , Tc. 
AU~ 0.3 mm/s; AV~ 50 pm/s 


- g-effects with crystal temperature variation, Tc = Tc(y). 
AU~ 0.3 mm/s; AV~ 43.75 pm/s 

Av velocity (m/s); (1 g 0 - 0 gj 

■ 9.045E-06 
R -8.312E-06 
R -4.303E-05 
R -7.774E-05 
R -1.125E-04 
R -1.472E-04 
R -1.819E-04 

■ -2.166E-04 


- g-sensitivity (horizontal case) based on max. buoyancy 
driven flow normal to growth direction is 10% of crystal 
growth rate (3mm/day or 0.035 pm/s) 
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transverse acceleration requirement: ~ lxl0‘ 4 g„ 




Gravity Effect on Velocity Difference (vertical 


- g-effects with crystal temperature variation, Tc = Tc(y). 
Stabilized orientation: AU~ 23.1 |am/s; AV~ 9.4 |j.m/s 

- g-effects with crystal temperature variation, Tc = Tc(y). 
Destabilized orientation: AU~ 18.1pm/s; AV~ 9.4 pm/s 

As far as transverse velocity difference is concerned both 
vertically stabilized and destabilized orientations have 
similar effects 

- g-sensitivity (vertical case) based on max. buoyancy 
driven flow normal to growth direction is 10% of crystal 
growth rate (3mm/day or 0.035 pm/s) 
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longitudinal acceleration requirement: ~ 2.7x10' g l} 




3-D Computational Grid and Code Validation 

Results 
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* Schiroky and Rosenberger (1984) 



ZnSe with Residual Gas - 3D Calculation Results 


• significant flow observed along the ampoule axis (z- 
direction) indicative of more deposition in the central area 
than near the walls 

• velocity contours in the cross planes(x-y) show appreciable 
variation only near the end walls (source and crystal) 

• species (Zn and Se 2 ) show fairly uniform distributions in the 
cross planes 

• predicted crystal growth rate from 2-D and 3-D calculations 
are in fair agreement 
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Summary of Theoretical Calculation 


• 2D and 3D calculations performed for ZnSe system 

• Residual gas effects considered 

• Calculations show that shear flow velocities of 10 to 50 
microns/s are induced by buoyancy effects (290 to 1400 
times growth rate) 

• g-level requirements established based on time scale 
analysis required transverse g level: < 1.2xl0' 4 g o 
required longitudinal g level < 8.5x10' g 0 

• It is noted that the Boussinesq model used in the 
calculations tend to underpredict velocities 
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